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Abstract. In this paper, steady-state MHD calculations 
of non-relativistic magnetized accretion discs driving jets 
are presented. For the first time, an energy equation de- 
scribing the effects of entropy generation along streamlines 
is included. Using a simplified approach, we showed that 
thermal effects have a tremendous influence on the disc 
capability to feed jets with mass. 

The disc ejection efficiency is measured by the parame- 
ter ^ = d\n Ma/d\nr, where Ma{r) is the local disc accre- 
tion rate. While previous self-similar solutions were only 
able to produce jets with ^ ~ O.OI, solutions with a coronal 
heating display huge efficiencies up to ^ '--^ 0.5. Moreover, 
taking thermal effects into account allows to obtain both 
fast and slow magnetic rotators. 

Since most of the jet properties (like asymptotic veloc- 
ity or degree of coUimation) depend on the mass load, it 
arises from this study that any quantitative result requires 
a detained analysis of the disc energetics. 

Key words: Accretion, accretion discs - Magnetohydro- 
dynamics (MHD) - Stars: formation - ISM: jets and out- 
flows - Galaxies: nuclei - Galaxies: jets 



1. Introduction 

1.1. Jets and discs 

Astrophysical jets are observed in a wide variety of ob- 
jects, from young protostars surrounded by a circumstellar 
accretion disc to compact objects in either, close binary 
systems or active galactic nuclei. In all these systems, sev- 
eral observational evidences point towards strong links be- 
tween accretion and ejection phenomena (Hartigan et al. 
1994 Mirabel et al. |1998| and Serjeant et al. [19981 ). Thus, 



theorists must provide models where accretion and ejec- 
tion can be simultaneously explained in a self-consistent 
way. Precisely, these models must provide a clear picture of 
how a fraction of the disc matter becomes an ejected flow, 
but also of the energetics (eg, the ratio of the jet kinetic 
power to the accretion power). Indeed, some problematic 



objects show powerful jets produced from discs showing a 
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much lower activity than expected (Rutten et al. 1992, Di 
Matteo et al. ^OOOD . We will call this the energy problem. 

It is well known that the "standard" accretion disc 
model of Shakura & Sunyaev ( 1973D cannot explain jets. 
In this model, all the mass is accreted towards the cen- 
tral object while the liberated mechanical power (called 
accretion power) is radiated at its surfaces. However, it 
as been recognized that, under certain circumstances, the 
disc could become geometrically thick, hence allowing for 
an efficient advection of the thermal energy (Abramowicz 
et al. 1988 ). This has an important observational conse- 
quence. Indeed, these advection-dominated accretion flow 
(ADAF) models allow a sizeable fraction of the accretion 
power to be actually advected onto the central object. 
Thus, such a power would be observationaly "lost" if this 
object happens to be a black hole. Moreover, it has been 
argued that such ADAFs would be hot enough to, hope- 
fully, provide a launching mechanism for jets (Narayan & 
Yi 1995). As a consequence, one might have jets with a 
high kinetic power arising from low-luminosity discs. All 
this reasoning is fine, but no comprehensive work has been 
done showing such a possibility. There is however an al- 
ternative to this energy problem. 



1.2. Models of magnetized discs driving jets 

First, we recall that a very efficient launching mecha- 
nism is the presence of a large scale bipolar magnetic field 
threading the accretion disc (for an alternative view, see 
Lery et al. 199£ ) . This magnetic field can be generated by a 
dynamo effect (see Rekowski et al. ^000 ) and/or by advec- 
tion of interstellar magnetic field (Mouschovias 1976). In- 
deed, since the seminal work of Blandford & Payne (1982), 
it has been showed that, combined with the disc rotation, 
such a field could (1) extract the disc angular momen- 
tum, (2) provide an acceleration force to a fraction of the 
mass and (3) maintain the jets self-collimated by a hoop- 
stress. Thus, explaining coUimated jets would require to 
revisit the standard model by self-consistently incorpo- 
rating magnetic fields. We call these magnetized objects, 
magnetized accretion-ejection structures (MAES). 
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The first attempt to describe such objects was made 
by Konigl (|1989| ) and followed by Wardle & Konigl ( |1993D , 
Li (1995, 1996). However, these works had two drawbacks. 
First, they did not self-consistently treat the disc verti- 
cal equilibrium, thereby obtaining huge parameter spaces 
and, second, the overall energy budget was not considered. 

A complete description of MAES has been undergone 
in a series of papers where, thanks to a self-similar ansatz, 
all dynamical terms were included in the magneto hydro- 
dynamic (MH D) e quations (Ferreira & Pelletier 1993a|, 



from the disc surface, the opening angle, defined as the 
angle between the poloidal magnetic field and the ver- 
tical axis, must be bigger than 30° (Blandford & Payne 



1995 , Ferreira |l997 , Casse & Ferreira 2000| ). In particular 



it has been shown that a necessary condition for jet pro- 
duction is the reversal of the sign of the magnetic torque 
at the disc surface: negative at the disc midplane (angu- 
lar momentum extraction), it must become positive at the 
surface in order to magnetically propell matter. This con- 
dition provides a strong constraint on the properties of the 
MHD turbulence believed to exist inside discs. Casse & 



Ferreira (2000) provided a comprehensive scan of the pa- 
rameter space of "cold" MAES, i.e. where enthalpy plays 
no role in jet dynamics. They assumed adiabatic magnetic 
surfaces and showed that, depending on the MHD turbu- 
lence properties, two kinds of "cold" MAES could exist: 

(1) discs where the magnetic torque due to the jets is much 
larger than the local "viscous" torque (of turbulent 
origin) . These MAES are called "non-dissipative" . In- 
deed, almost all the available power feeds the jets and 
only a tiny fraction is locally dissipated and radiated 
by the disc (see also Ferreira & Pelletier 1993a and 
Ferreira 1997). Another important consequence of the 
larger magnetic torque is a disc density smaller than in 
a standard disc (for the same accretion rate). A non- 
dissipative MAES could then be optically thin, and is 
thus expected to provide a spectrum very much com- 
parable to that of an ADAF. 

(2) discs where the magnetic torque is comparable to the 
viscous torque, called "dissipative" MAES. Such struc- 
tures dissipate as radiation a power comparable to that 
carried away by the jets. No solution has been found 
with a dominant viscous torque. 

Thus, MAES quite straightforwardly address the en- 
ergy problem mentionned earlier. Instead of advecting en- 
ergy onto the central object like ADAF do, MAES put a 
sizeable fraction of the accretion energy into jets. In fact, 
everything depends on the properties of the MHD turbu- 
lence inside the disc, which are still unknown. 

1.3. The "cold jet" approximation 

Nonetheless, all these results (as well as all previous works 
on disc winds) rely on what happens to be a very impor- 
tant assumption, namely the "cold jet " approximation. 
Usually, in the literature, a "cold jet" is a jet where en- 
thalpy plays no role in its launching (ie, is negligible in 
the Bernoulli equation). Thus, in order to get ejection 



1982). However, the crucial question is the amount of mass 
that can be loaded in the jet. This is never addressed in 
any steady-state numerical simulation, since the underly- 
ing accretion disc is merely a boundary condition. This 
mass load critically depends on the force balance at the 
disc surface, where matter can still cross the magnetic 
field lines. In this transition layer from a resistive disc to 
an ideal MHD jet, the only force allowing matter to go 
up is the plasma pressure gradient (Ferreira & Pelletier 
1995). Thus, the vertical profile of the temperature does 



play a role in the launching process. Since it mainly occurs 
in the resistive MHD region, a negligible enthalpy in the 
ideal MHD zone does not mean that it had no influence 
on jet formation. 

In this paper, we will show that the precise treatment 
of the energy equation must be done correctly, especially 
at the disc surface. This is a sine qua non condition in 
order to get realistic quantitative answers concerning the 
capability of accretion discs to launch jets. In particular, 
the presence of a hot corona could drastically change the 
whole picture. Because cooling and heating processes de- 
pend on density and temperature, the conditions at the 
disc surface vary with the radius and depend also on the 
astrophysical context studied (e.g around a compact ob- 
ject or a young star). Here, our purpose is just to illustrate 
this stinging but crucial point by using some crude pre- 
scription. We will investigate specific astrophysical con- 
texts in forthcoming papers. 

In section 2, we recall the basic stationary MHD equa- 
tions describing a MAES and define the relevant disc and 
jet parameters. Section 3 is the cornerstone of this pa- 
per. We introduce the local disc energy equation with the 
prescription we used and discuss the constraints imposed 
by the global energy conservation. Thermal effects on sta- 
tionary jet production are then analyzed. In section 4, we 
show some typical self-similar MAES solutions and display 
two extreme cases. We finally conclude in section 5. 

2. Magnetized accretion-ejection structures 

2.1. Basic MHD stationary equations 

Thanks to axisymmetry, the vectorial quantities expressed 
in cylindrical coordinates (r, (/), z), can be decomposed into 
poloidal and toroidal components, e.g. u — Up + flre^ 
and B = Bp -t- B^Bcf,. A bipolar configuration allows us 
to describe the poloidal magnetic field as 



Bj, — -Va X Ed 

r 



(1) 



where a(r, z) is an even function of z and a = constant 
describes a surface of constant magnetic flux (a — rA^, 
Aff, being the toroidal component of the potential vector). 
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The following set of equations describes a steady-state, 
non-relativistic MAES, 

— Mass conservation 

V • pit = (2) 

— Momentum conservation 

pu-Vu^ -VP - pV$G + J X B + V • T (3) 

— Ohm's law and toroidal field induction 

TlmJtpe^ = UpX Bp (4) 

V • i'^VrB^) = V • -{B^Up - Bpflr) (5) 



where p is the density of matter, P the thermal pres- 
sure, $G — —GM^,/{r^+z^y/^ the gravitational potential, 
J = V X B / po the current, T the "viscous" stress ten- 
sor, — rjni/po and v!^ the (anomalous) poloidal and 
toroidal magnetic diffusivities. All transport coefficients 
appearing in the above set of equations, namely Vm, i^'m 
and ( "viscosity" , contained in T) , are assumed to be 
of turbulent origin. The state equation used is the perfect 
gas law 



P = nkuT 



(6) 



where A;b is the Boltzmann constant, n = p/rup (nip being 
the proton mass) and T the plasma temperature. 

The last equation is the energy equation. Until now, it 
has been replaced by a simple prescription on the temper- 
ature, namely isothermal or adiabatic magnetic surfaces 
(see eg Ferreira & Pelletier 1995 for a discussion on this 
subject). The aim of the present paper is to study the ther- 
mal effects on jet production. Therefore, we must consider 
the energy equation, which can be written as 



pT^ = pTup ■VS = Q 
at 



(7) 



where S is the specific entropy and Q is the local source of 
entropy, arising from the difference between heating and 
cooling processes (see Sect. 3 for more details). 

2.2. Dynamical parameters of MAES 

The set of equations (||) to (^ completely describes a 
steady-state, non-relativistic MAES. We recall here the 
definitions of the parameters used: 



The disc aspect ratio 



(8) 



The MHD turbulence level (measuring the strength of 
the poloidal magnetic diffusivity at the disc midplane) 



(9) 



— The magnetic Prandtl number (measuring the strength 
of the turbulent viscosity) 



-p - — 

Vrr,. 



z = 



The magnetic diffusivity anisotropy 



(10) 



(11) 



where h{r) is the vertical disc scale height and Va the 
Alfven speed at the disc midplane. These four parameters 
are defined locally and could, in principle, vary along the 
disc. In our study however, they are assumed constant. 
Also, while the disc aspect ratio e remains free (because 
we did not treat the real energy equation) , the three other 
turbulence parameters must remain free, until the devel- 
opment of a full theory of MHD turbulence. All other disc 
and jet quantities can be expressed in terms of these pa- 
rameters (see Casse & Ferreira 2000 , hereafter CF). As an 
example, the usual alpha coefficient of Shakura & Sunyaev 
is given by 



Vy 



V n* 



(12) 



where a*^ = v^/Csh = a,„(/i/7)i/2, p = {VA/^Kh)"^ is a 
measure of the magnetic field strength (JIk is the Keple- 
rian rotation rate) and 7 the adiabatic index. Note that, 
since steady state MAES require fields c lose t o cquipar- 
tition (ie. p < 1, Ferreira & Pelletier 1995 ), one has 



2.3. Generalized ideal MHD integrals 

Inside the disc, accreting matter must cross the field lines. 
Thus, a stationary configuration requires turbulent mag- 
netic diffusivities. This turbulence is believed to be trig- 
gered inside the disc by some MHD instability, but decays 
vertically on (say) a disc scale height. In the ideal MHD 
regime characterizing the jets, integrals of motion appear 
in the basic MHD equations. These are 



(13) 
(14) 
(15) 



Up 


via) 


Bp 




l^oP 


ri*(a) 








77 




Popr 




= nr^ 


rB^ 




77 



where 77(a) = y^PoPA is a constant along a magnetic sur- 
face and PA is the density at the Alfven point, where the 
poloidal velocity Up reaches the poloidal Alfven velocity 
VA.p. The rotation rate of a magnetic surface is ri,(a) 
and VA is the Alfven radius. Following Blandford & Payne 
(1982), we introduce the following jet parameters 



z=0 



A = 



Q r-2 
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^qTq 

Vap,a 



(16) 



The subscript "o" refers to quantities evaluated at the disc 
midplane, "A" at the Alfven surface. These parameters 
are directly related to the disc parameters (see CF). The 
last ideal MHD invariant is the Bernoulli integral. This in- 
variant is modified by the presence of heating, contrary to 
the three former invariants. By projecting the momentum 
conservation on a magnetic surface and usin g ther mody- 
namics relations, one gets (Chan & Henriksen 198C| , Sauty 
& Tsinganos 1994) 

E{a) = — + —^- + 'PG~n,^-nis,a) (17) 

where the enthalpy (second term) will be hereafter repre- 
sented by H and 



H(s, a) 



Q(s',a) 



+ p{s' , a)up{s' , a) 



ds' 



(18) 



is a heating term that depends on the curvilinear co- 
ordinate s along a given magnetic surface (labelled by 
a constant a{r,z)). Here, s+ represents the coordinate 
of the transition between the resistive and the ideal 
MHD regimes (above the disc surface). Jets from quasi- 
Keplerian accretion discs are energetically possible if 
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7-1 To 



(19) 



is positive at the base of the jet. Thus, either the magnetic 
lever arm A is larger than 3/2 or some heating occurred in 
the underlying resistive layers. This last possibility arises 
whenever the disc is already hot (ie. thick e ~ 1 or slim 
e < 0.3 discs) and/or if some coronal heating leads to an 
increase of the temperature T"*" at the disc surface. 

3. Entropy generation in MAES 

3.1. Local entropy generation 
The local energy equation ^ is 



Q = pT 



dS 



dt 

7 ki 



dH 
^~dt 



dP 
'dt 



7 - 



1 m. 



■pUp ■ VT - Up ■ VP 



(20) 



In a thin (or slim) accretion disc, the gradient of any 
plasma quantity is mostly vertical. So the main influence 
of an entropy source Q occurs whenever the plasma veloc- 
ity becomes vertical namely, at the disc surface. There, a 
non- vanishing Q will have two distinct and important con- 
sequences. First, it raises the plasma temperature (T+) 
and, second, it causes an increase of the vertical plasma 



pressure gradient. Now, the vertical component of the mo- 
mentum conservation writes 



puz 



duz 
dz 



dP 

dz 



— pfl%z — 



dBl 



Bl 



dz 2p,o 



(21) 



In this equation, the only force that can gently expell mat- 
ter from the disc is the vertical plasma pressure gradient. 
Thus, increasing this gradient through an entropy source 
Q allows new configurations to exist. For example, with 
parameters (e, am, Vm, Xm) similar to those of a "cold" 
MAES (ie. Q — 0), a solution with an entropy source 
would display a larger vertical force, allowing thereby a 
larger mass flux k. Or, with the same entropy source, one 
could get a steady ejection from a MAES with a larger 
magnetic pinching, corresponding to turbulence parame- 
ters impossible to achieve with Q ~ 0. 

Thus, we see that both the amplitude and the vertical 
profile of the entropy source Q are important features for 
stationary MAES. In fact, there is a local entropy source 
Q only if there is a slight discrepancy between the local 
heating and cooling terms, namely 



Q — (fe// + I'turb + Fext) " {J^rad + Murb) 

where 



(22) 



(23) 



is the effective Joule and viscous heating, Tturt is a tur- 
bulent heating term which cannot be described by simple 
anomalous transport coefhcients (eg, the term Teff), T^xt 
is an external source of energy (typically due to some il- 
lumination by UV or X-rays, hence maximum at the disc 
surface) , 



Arad — V • 5'. 



rad 



(24) 



is the radiative cooling {Srad being the radiative flux) and 
Aturb is a cooling due to a turbulent energy transport, 
which is most probably also taking place inside turbulent 
discs (see eg. Shakura, Sunyaev & Zilitinkevitch 1978). 

If the disc is optically thick, the radiative flux Srad 
can be written in a simple way using the diffusion approx- 
imation (eg. Ferreira & Pelletier 1993a). However, both 
the amplitude and the vertical profile of the turbulent en- 
ergy transport Aturt arc unknown. This is really problem- 
atic, since this term could be the leading cooling term 
(like convection, see Ferreira & Pelletier 1995). On the 
other hand, magnetized discs may well have active coro- 
nae, hence a non negligible Tturb mainly acting at the 



disc surface (Galeev et al. 1979, Heyvaerts & Priest 1989, 
Miller & Stone 2000). Therefore, even if we neglect any 



external influence, the only term whose amplitude and 
vertical profile are re ally known is T^ff. Thus, following 
Narayan & Yi ( 1995 ) for ADAFs, we define a parameter 
/ as 



UQdV JyQdV 



Iv^effdV P. 



diss 



(25) 
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where the integral is made over the volume V occupied 
by the disc and Pdiss is the total power which could be 
dissipated inside the disc. With this prescription, the am- 
plitude of the entropy source Q is controlled by /, but not 
its vertical profile. This profile depends on the relative im- 
portance of the various terms appearing in Eq. (|2|) and, 
as such, is beyond our actual knowledge. To tackle this 
difficult problem, we will assume the profile and discuss 
the corresponding physical implications on jet production. 

3.2. Global energy conservation 

A strong constraint on this parameter / arises from the 
global energy conservation. If we integrate the local energy 
conservation equation over the whole volume V, we get the 
following energy budget 



Pn. 



(26) 



where Pace is the power liberated through the accretion 
process. 



Pext — / ^extdV 
'V 



(27) 



is the total input of energy due to an external source, 
Pjet is the total power (mechanical, magnetic and ther- 
mal) carried away by one jet (hence leaving the disc at its 
surface) and Prad the luminosity radiated at one surface 
(see Appendix A). 

In this paper, we do not consider any external source 
of energy, thus Pext = 0. Moreover, we do not expect in a 
quasi-Keplerian flow the accreted turbulent energy to be 
significant with respect to the mechanical energy. Thus, 
we assume that the total turbulent power is negligible, 
namely 



Pturb = / {^turb — ^turb)dV — 

Jv 



(28) 



The characteristic power, liberated by a MAES located 
around a central mass , between the radii and r^, and 
fed with an accretion rate Mae = Afa(re), is 



Piib = rjii 



GMMa 

' 2^ 



(29) 



where the efficiency rjub is given in Appendix B. In discs 
producing jets, the accretion rate is a function of the ra- 
dius. We therefore define the local ejection efficiency as 



d In Ma 
dlnr 



(30) 



In terms of the fiducial value Pm, , the global energy con- 
servation (^6|) of thin (or even slim) MAES writes (see 
Appendix B for details) 



Pace , / 1 eA 

"'^'^ (1 - ^) 1 + 



Pih 



2 1 + A 



(31) 



providing jets powered with 



2P, 



jet 



A 



lib 



1 + A 



27 T+^^, 



and a disc luminosity of 



IP 



p 



lib 



(1-/) 



Pdis 

~pr 



Pn. 



1 T„ 



2P. 



jet 



lib 



P 



lib 



(32) 



(33) 



Here, the parameter A is the ratio of the magnetic torque 
to the viscous torque at the disc midplane and qis a mea- 
sure of the radial current density at the disc midplane (see 
CF for more details). 

From the above energy budget, it becomes obvious that 
two extreme cases can be accounted for, namely 

— / <C 1: cooling and heating processes balance each 
other very effectively at every point, so that there is 
no entropy generation at all (Q — 0). In this case, the 
dissipated power Pdiss is finally radiated at the disc 
surfaces. 

— / = 1: a large entropy source is generated in this case, 
because of a turbulent transport more efficient than 
radiative cooling and the presence of another heat- 
ing process (eg. coronal heating). As a consequence, 
no power is radiated at all, the accretion power being 
transferred into the jets. 

The case / > 1 can still be obtained but requires an extra 
source of energy {Pext or Pturb)- This is also true when 
/ < 1 but the disc luminosity exceeds that provided by 
Eq. (H. 

3.3. Constraints on steady-state jet production 

The trans-Alfvenic condition is modified by the possible 
presence of a heating term along the flow. This can be 
directly seen in the function 5, deflned as 17 = r2*(l — g) 
and evaluated at the Alfven surface 

3 1 _^ 2 
2" ~^ 



9\ 



= 1- 



A Lo\ A3/2(l + 4/r^)i/2 
His A, a) 



7 



T+ 



7- 



1 







(34) 



The subscripts "o" and "A" denote a quantity measured, 
respectively, at the disc midplane and at the Alfven sur- 
face. A necessary condition for trans-Alfvenic cold jets is 
usually expressed as 



LJ\ > 1 



(35) 



implying t hese je ts are fas t rota tors (Michel 1969 , PcUetier 
& Pudritz 1992, Ferreira 1997 ). However, this condition 
is only strictly valid for very large magnetic lever arms 
(A — !■ 00). The generalized trans-Alfvenic condition is in 
fact 



9a >0 



(36) 
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Thus, Eq. (|3J) shows that entropy generation inside the 
disc (leading to a large T+) or alternatively, some coronal 
heating (significant 7i(s^, a)), allows the fulfillment of the 
above condition with to a < 1- Therefore, we expect slow 
rotators to emerge as (magneto-)thermally driven jets. 

Following the same reasoning than Ferreira (1997) in 
Sect. 3, but considering the generalized Bernoulli equa- 
tion, we obtain an inequality constraining the magnetic 
lever arm 



A 



A- 1 



< A 



Al/2(1 + 4/^2^)1/2 



7- 



1 Tn 



(37) 



This quite complex relation is equivalent to A > Xmin — 
1 + C- This shows that, to produce super- Alfvenic jets, 
there must be a minimum magnetic lever arm Xmin- 
Namely, steady-state MAES cannot display an ejection 
efficiency larger than 

Cma. = 2CTTA ^^^^ 

The exact value of C is strongly dependent on Q (both 
amplitude / and vertical profile) . Basically, the more heat 
is added and the more mass is ejected. For example, for 
T+ <To, C = 5 while for T+ ~ s-^To, C ^ 0.3. Thus, as 
already pointed out, the major constraints on ^ arise from 
the disc vertical equilibrium in the resistive MHD regime 
and the energy balance. 



4. Numerical results 

The self-similar ansatz 

In this section, we briefly describe the radial self-similarity 
that enables us to consider both the resistive, viscous disc 
and the ideal MHD jet, including an exact treatment of all 
dynamical terms. Any physical quantity A{r, z) is written 
as 



- ) Fa{x) 



(39) 



where x = z/h(r) and is the outer radius of the mag- 
netized disc. This special ansatz allows to separate the 
variables so that the full MHD equations become a set of 
non-ordinary (they contain three singularities) differen- 
tial equations. The values of the radial exponents (given 
in Ferreira & Pelletier 1993a ) arc all linked to the ejection 
index ^. For example, the magnetic flux becomes 

a{r, z) ^ aJ -Y Fa{x) (40) 



with p = 3/4-1-^/2. Inside this framework, the entropy 
source Q must be written as 

Q^qJ^\ Fgix) (41) 




x-z/h 



Fig. 1. Vertical profile of the entropy generation function 
Q for a given ejection index ^ — 0.1. The value of this 
function at the disc midplane {x — 0) is imposed by the 
self-similar radial scalings (see text). The integral of this 
function is measured by the parameter /. The particular 
shape shown here (increasing towards the disc surface) 
mimics a coronal heating. 



where, using Eq. (20), its value at the disc midplane must 
satisfy 



Qo = Q(r,z = 0) = 



7 



7-1 
PqUpq 



ap 



PoU. 



o ^po 



(42) 



This is a very strong constraint which also occurs in self- 
similar ADAF models. There ^ = 0, so an ADAF must 
have 7 7^ 5/3 in order to get some advected heating at the 
disc midplane. Here, we choose an adiabatic index 7 = 5/3 
and scale the vertical profile Fq{x) with the parameter 
/. We display in Fig .(||) a typical profile corresponding 
to a discrepancy between cooling and heating processes, 
occuring mainly at the disc surface. We did however play 
with slightly different vertical profiles (see Fig. (^). 

The set of MHD equations has three critical points 
in the ideal MHD regime. These critical points are the 
slow-magnetosonic, the Alfven and the fast-magnetosonic 
points. The solutions presented here share the same be- 
havior than in previous studies (Ferreira 1997 , CF). The 
magnetized flow passes through the slow and Alfven points 
(the regularity conditions impose the values of /x and ^) 
and then encounters the fast-magnetosonic point in a rec- 
oUimation motion. None of these solutions pass through 
the last critical point (see Vlahakis et al, in preparation). 

4-.2. From magnetically-driven to thermally -driven jets 

As showed in Sect 3.1, the entropy generation increases the 
plasma pressure gradient, thereby allowing an enhanced 
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f f 

Fig. 2. Variation of several quantities with the entropy 
parameter / for a given magnetic configuration: e = 0.1, 
am = 1-5 (a*j — 0.8), Vm = 1, Xm = 1-5. Non-relativistic 
jets from Keplerian accretion discs can be either fast ro- 
tators [ujA > 1) or slow rotators {uja < 1), depending 
on the amplitude of /. The coronal temperature increases 
monotonously with /. The plots in mass load k and mag- 
netic lever arm A illustrate the fact that the more entropy 
is generated the more mass is loaded (see Sect. 3.1). 



vertical mass flux. In particular, jets are now possible in 
configurations where the magnetic pinching force would be 
overwhelming and squeeze the disc. Thus, thermal effects 
significantly enlarge the parameter space of MAES (e, a™, 

'Pmj Xm)- 

In order to illustrate such an effect (strong magnetic 
compression balanced by entropy generation), we choose 
a MAES configuration described by the following param- 
eters: e = 0.1, am = 1-5 {a'^ = 0.8), Vm = 1, Xm = 1-5. 
The chosen disc aspect ratio is relevant in the astrophys- 
ical systems we consider (see Appendix C). Figure (||) 
shows relevant quantities obtained with the same verti- 
cal profile Fq (see Fig. (^) and increasing values for /. 

First, a minimum heating (/ ^ 0.001) is here required 
otherwise the thermal plasma gradient cannot counter- 
act the magnetic pinching force. Note that this mini- 
mum heating decreases with decreasing am (because of 
the weaker toroidal magnetic compression). Second, as ex- 
pected, both the mass flux k and coronal temperature T'^ 
grow as / increases. Their maximum values depend here 
on the chosen set of disc parameters. In a consistent way, 
the magnetic lever arm A diminishes with /. Third, as / 
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Fig. 3. Bernoulli invariant E{a) and its components along 
a given magnetic surface (s(a) = z{a)lho)-, for the hottest 
jet (/ = 0.82) in Fig. (||). The quantities are normalized to 
the square of the angular velocity at the footpoint of the 
magnetic surface (il^r^). Enthalpy is of the same order 
than the MHD Foynting flux in the corona, signature of a 
"hot" jet. 
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Fig. 4. Same plot than in Fig. (^), but for the coldest jet 
(/ = 0.002) of Fig. d). The enthalpy does not play any 
role in the jet energetics. This is the signature of a "cold" 
jet. 

increases, thermal effects become significant and there is 
a transition between "cold" and "hot" outflows. Indeed, 
when / is very small, the jets show the usual character- 
istics of a "cold" outflow, namely a fast magnetic rotator 
(iOA > 1), a large magnetic lever arm A, a small mass load 
K and a coronal temperature smaller than that in the 
disc. At the opposite, when / is close to unity, the jets 
are slow magnetic rotators (w^ < l)i have small magnetic 
lever arms, high mass loads and high coronal tempera- 
tures. 

We will label as "hot" , jets whose initial enthalpy plays 
a significant role in their launching. We display in Figures 
(H) and (^), the Bernoulli constants E{a) and their dif- 
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Fig. 5. Ratio I a/ Ism of the current at the Alfven point 
to the current at the slow-magnetosonic point as a func- 
tion of the entropy parameter /, for the same MAES 
as in Fig. (||). As / increases, there is a transition from 
one regime to the other. The first regime corresponds to 
"cold" jets, where the thermal energy is insignificant and 
the mass load is small. The second regime corresponds 
to "hot" jets, where a large thermal energy supplants the 
Poynting flux in propelling matter in the sub-Alfvcnic re- 
gion. 



ferent components, for the two extreme jets obtained in 
Fig. (|) (respectively / = 0.82 and / = 0.002). In the 
"hot" jet, the enthalpy at the base is of the same order 
than the MHD Poynting flux. The jet acceleration in the 
sub-Alfvenic region is mostly due to the decrease of the 
enthalpy. This shows the importance of thermal effects, 
although the Blandford & Payne criterion (opening angle 
larger than 30°) is, here, still verified. In the "cold" jet 
however, the enthalpy is negligeable and the acceleration 
is due to the decrease of the poloidal current. 

In order to clearly distinguish "cold" (magnetically- 
driven) from "hot" (thermally-driven) jets, we look at the 
decrease along the jet of the poloidal current, defined as 



27r 

/ = —rB4 

fJ.o 



(43) 



Indeed, this current is a measure of the transfert of MHD 



Poynting flux into kinetic energy flux (see Ferreira 1997). 
We plot in Fig. (||) the ratio I a/ Ism of the current still 
available at the Alfven point to the current provided at 
the base of the flow (at the slow-magnetosonic point). For 
small / (here / < O.I), this ratio remains high because 
a tiny fraction of mass has been loaded in the jet. As / 
increases, this ratio decreases because the enhanced mass 
load requires more magnetic energy to be used. But at 
some point (/ > 0.1), the thermal energy reservoir sup- 
plants the MHD Poynting flux in propelling matter until 



the Alfven point. Thus, the ratio is now increasing be- 
cause of this new energy supply. The two regimes shown 
in Fig. (||) define a "cold jet" and a "hot jet" zone. 
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Fig. 6. Same plots than in Fig. (g), but for a weaker turbu- 
lence level am = 1.2. The two jet regimes are now splitted. 
"Cold" jets exist only for / < 0.02, "hot" jets only for a 
large entropy generation, / > 0.5. 

For a given MAES, there is always an upper limit to 
the entropy generation (besides the overall energy conser- 
vation constraint) coming from the requirement that gA 
must always remain smaller than unity (ie. remains 
negative). Indeed, the available current at the Alfven point 
writes 



Ia 
Ism 



X 



gA 



A- I 



(44) 



The increase of this ratio for increasing / is mainly due to 
the increase of qa (see Fig. (||)). Thus, there is a maximum 
/ allowed (which can be smaller than unity) for steady- 
state MAES. 

A transition from "cold" to "hot" jets requires to cross 
a zone where a large mass load (provided by the enhanced 
plasma pressure gradient) is still accelerated. This is only 
achieved if a sufficient amount of toroidal magnetic field is 
present {B^ is the most relevant ingredient of the Poynt- 
ing flux) and/or if the thermal energy reservoir is enough. 
However, depending on the magnetic configuration, the 
jet mass load could increase with / faster than the en- 
thalpy. As a result, configurations with a weak toroidal 
magnetic field would not be able to steadily launch jets in 
a zone where the enthalpy is not sufficient. As an exam- 
ple, we provide in Figures (^) and (^) a magnetic con- 
figuration with a smaller turbulence level (a™ = 1-2), 
hence a toroidal magnetic field smaller than in the pre- 
vious case. This time, no solution has been found between 
0.02 < / < 0.5. 
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Fig. 7. Same plot than in Fig. (||), but for a MAES with 
a weaker turbulence level am = 1-2. A forbidden zone ap- 
pears where both the MHD Poynting flux and the thermal 
energy are too small to propel the mass loaded. 

4.3. Two extreme MAES 

Because of the extreme sensibility of the mass load k to the 
entropy source Q, we will not give a parameter space in- 
cluding these thermal effects. Indeed, both the amplitude 
/ and the vertical profile Fq should be obtained by the 
resolution of a realistic energy equation including cooling 
and heating terms. These terms (especially the cooling 
function) depend critically on the astrophysical context 
studied. Our purpose here is to show the broad range of 
available solutions by providing two extreme MAES con- 
figurations. 

4.3.1. A dense, slow "hot" jet 

In this extreme solution, a huge (/ = 1) entropy genera- 
tion occurs above the disc surface, see pannel e) in Fig. (^. 
This might be achieved with some sort of non-local coro- 
nal heating, efficiently transporting upwards the thermal 
energy released in the underlying layers of the disc. 

For this pedagogical example, we choose a magnetic 
configuration described by the following disc parameters: 

£ = 0.1, a,„ = 4 «„ ^ 1.8), Vrn = 1, - 1.4. It 

exhibits straight poloidal magnetic field lines inside the 
disc and a factor q = 1.012, leading to a torque ratio 
A — 2.6 (comparable magnetic and viscous torques). The 
required magnetic field must be quite small, fi = 0.26 (ie. 
a plasma beta of 7.7 at the disc midplane). 

The enormous release of thermal energy allowed a fan- 
tastic ejection, namely ^ = 0.456. The corresponding jet 
parameters are a mass load k = 1.04, a magnetic lever 
arm A — 1.9 and a rotator parameter lua — 1.65. Al- 
though this solution is "hot" , it is still a fast rotator. This 
is because the toroidal field at the disc surface is high. 
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Fig. 8. Plots of various characteristics of the two extreme 
("hot" and "cold") jets presented in Sect. 4.3: a poloidal 
cross-section of the magnetic surfaces anchored at a radius 
To; b poloidal velocity of the plasma along a magnetic sur- 
face in units of the Keplerian speed at its footpoint; c and 
d, plasma density and temperature along a magnetic sur- 
face normalized to their values at the disc midplane. The 
"hot" , dense jet is drawn in solid line while the "cold" , 
tenuous jet is in dashed line. The cross symbolizes the lo- 
cus of the slow-magnetosonic (SM) point and the star the 
locus of the Alfvcn (A) point. Pannels e and f show the 
effective (viscous and Ohmic) heating term Teff and the 
prescribed entropy source Q at a constant radius, normal- 
ized to the same fiducial quantity. 

I — 0.94i?o, despite straight poloidal field lines (large 
a„i, see CF). 

The disc corona is hot and quite dense and the jet does 
not widen much before recoUimation occurs (see Fig. (H)). 
This widening is however sufficient to produce jets where 
most of the available energy has been transferred into 
poloidal kinetic energy, namely 



up^ ~ ^2E{a) 



(45) 
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along every magnetic surface. The presence of a large 
thermal energy allows a terminal velocity larger than 
r^o'^o •\/2A — 3, which is the characteristic value obtained 
in "cold" flows {flora is roughly the Keplerian speed at 
the footpoint ro of the magnetic surface). Figure (^) shows 
a cross section of such a solution, settled around a com- 
pact object (upper image). Even in this extremely "hot" 
solution, the jet opening angle at the transition between 
resistive and ideal MHD regimes is 6^ = 37", thus still 
fulfilling the Blandford & Payne condition for "cold" jets. 

The global energy budget is here exceedingly simple. 
Since / = 1, all the accretion power is finally converted 
into jet kinetic power. Pace = 2Pjet — 0.58 Put- Again, we 
do not claim this is a realistic solution, but we want to 
point out at the extreme sensibility of the produced jets 
to the disc energetics. 



4.3.2. A tenuous, fast "cold" jet 

Here, we show an example of a jet launched from a MAES 
with an almost negligible entropy generation, namely / = 
5 lO^'^ see pannel f) in Fig. (^). This slight discrepancy be- 
tween heating and cooling processes occurs only inside the 
disc. Thus, this solutions requires no energy source above 
the disc, neither a coronal heating nor an illumination. 

The MAES is described by the following disc parame- 
ters: £ = 0.1, am = 0.3 « = 0.2), Vm = 1, Xm = 0.047. 
It exhibits a very strong curvature of the poloidal magnetic 
field lines inside the disc and a factor q = 0.25, leading to 
a torque ratio A — 15.67 (dominant magnetic torque at 
the disc midplane). Here, the required magnetic field is 
closer to equipartition, /i = 0.745 (ie. a plasma beta of 2.7 
at the disc midplane). 

In this case, the strong compression due to the vertical 
magnetic pressure gradient almost forbids the formation of 
a jet. The presence of a slight entropy generation however 
allows an ejection with an efficiency ^ — 0.001. The corre- 
sponding jet parameters are a tiny mass load k — 2.8 lO"**, 
a huge magnetic lever arm A = 440 and a rotator param- 
eter u!A — 1.2. Again, this is not a "very" fast rotator 
despite the small mass load, because of the small value of 
the toroidal field at the disc surface {\B'^\ = 0.12_Bo). 

The jet is "cold" , with a very strong adiabatic cooling 
due to its large widening (see Fig. (^). This widening re- 
sulted in a very efficient acceleration of the ejected plasma, 
leading to an asymptotic poloidal velocity ~ 23 times the 
Keplerian speed at the field lines footpoint. This is because 
of the huge magnetic lever arm A. In particular, relativis- 
tic speeds could be expected if such a MAES is settled at 
the vicinity of a compact object. Indeed, if we define CToo 
as the asymptotic ratio of the MHD Poynting fiux to the 
kinetic energy flux, one gets an asymptotic Lorentz factor 
of 



where rg = 2GAhJc^ is the Schwarzschild radius of the 
central object. Of course, a full calculation of the relativis- 
tic MHD equations should be undergone in order to get 
the value of (Too- As an example, if some equipartition is 
achieved (aoo = 1, Li et al. 1992D , one could in princi- 
ple get jets with asymptotic Lorentz factors ranging from 
~ 36 [to ~ 3rg) to - 4 {r^ ~ 30rg). 

Figure (^ shows a cross section of such a solution, 
settled around a compact object (lower image). The open- 
ing angle 6'+ is quite large, 53°, typical of a "cold" out- 
flow. The global energy budget is here more subtle than in 
the previous case. Indeed, we obtain an accretion power 
of Pace — 1-09 Piib, jets fed with a total power of only 
2Pjct — 0.45 Piib with the remaining power being radi- 
ated at the disc surfaces, namely 2Prad — 0.64 Puf,. Thus, 
we do have a dissipative MAES producing "cold" jets. 
The very strong curvature of the poloidal magnetic fleld 
lines yields a strong toroidal current, hence a large Ohmic 
dissipation. Moreover, the large anisotropy of the MHD 
turbulence {xm = 0.047), provides also a large dissipation 
of the toroidal magnetic field (see Eq. (p3|)). Hence, the 
dissipated power Pdiss, which is here converted into radi- 
ation, is quite large. Equation ( ^^ shows that the exact 
amount of energy powering the jets depends on the ra- 
tio \B'^/qBo\. Here this ratio is equal to 0.48, namely a 
generated toroidal magnetic field at the disc surface half 
the expected value (measured by q). This decrease of the 
toroidal magnetic field is due to this enhanced magnetic 
diffusivity v'^. 

5. Conclusion 

In this paper, we removed a basic assumption of most 
theoretical works dealing with axisymmetric, stationary 
Keplerian discs launching jets (MAES). Indeed, the disc 
energy equation was usually neglected, using either adi- 
abatic (Casse & Ferreira 2000) or isothermal (Wardle & 



1 



2A- 3 

12(1 + (Too) 



3^0 



(46) 



Konigl |1993| Li |1995| |1996| , Ferreira |1997| ) magnetic sur- 
faces. 

We used a self-similar ansatz allowing us to take into 
account all dynamical terms in the resistive, viscous disc 
as well as in the ideal MHD jet. The transition between 
the accretion and the ejection motions is self-consistently 
described. We modelled the entropy increase along each 
streamline by prescribing the vertical proflle of the heat 
effectively converted into entropy. This is an approach very 
similar to that used in ADAF models, with a signiflcant 
distinction. Indeed, in the case of MAES, the main part 
of the heat released in the disc could be absorbed by the 
jet instead of being advected onto the central object. 

The possibility to steadily produce jets lies in the 
rather subtle disc vertical equilibrium. Both gravity and 
the chosen bipolar magnetic topology vertically pinch the 
accretion disc. Thus, mass can only be lifted from the 
disc by the plasma pressure gradient. Because a local en- 
tropy generation enhances this gradient (see Eq. (^o|)). 
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Fig. 9. Cross-sections of two extrenie magnetized accretion-ejection structures (MAES), settled around a neutron star 
of Af, = IAAIq and fed with an accretion rate of Mae — 10^^^ MQ/yr. The poloidal magnetic field is shown in solid 
lines, streamlines are dashed and the grey contours represent the logarithm of the plasma temperature. Both MAES 
have the same temperature at the disc midplane. The upper image is an extremely "hot" MAES (/ = 1) with an 
entropy source located above the disc surface. This solution has a huge ejection efficiency of ^ = 0.456, plasma reaching 
an asymptotic velocity 1.6 times the Keplerian speed flora at the magnetic field lines footpoint Tq. The lower image is 
a "cold" MAES (/ = 5 10~^) with no heating source above the disc. Such a solution has ^ = 0.001 with asymptotic 
velocities 23 times fio'^o- 



the amount of mass escaping the disc can be much larger. 
Alternatively, tiny mass loads, impossible to achieve in 
previous MAES calculations, are now reachable. Indeed, a 
huge magnetic compression can be now balanced by some 
amount of heating. Moreover, there is a second important 
consequence of entropy generation, namely enthalpy gen- 
eration. This enthalpy is carried away by the outflowing 
material and can be converted afterwards into kinetic en- 
ergy. To summarize, entropy generation results in dynam- 
ical (measured by the ejection index ^) and thermal (mea- 
sured by the coronal temperature T+) effects. For the first 
time, this is self-consistently illustrated in a global energy 
budget of MAES. This paper has illustrated that heating 
is as important in steady jet formation from a accretion 
disc as in steady outfl ow release around embedded sources 
(e.g. Lery et al. 



19991) 



We are now in a position to provide clear diagnos- 
tics on steady-state MHD simulations of jets from Keple- 
rian accretion discs. Indeed, all effects acting on the ejec- 
tion process itself are taken into account. Of course, this 



study was done under the self-similar assumption, but the 
links between MHD invariants are local and, as such, in- 
dependent of this ansatz. Besides, in a quasi-Keplerian 
disc where gravity (which is self-similar) is the domi- 
nant driving term, we expect self-similar solutions to be 
quite reasonable. For example, local jet parameters (mass 
load K, magnetic lever arm A and rotator parameter coa) 
as well as the shape of the Alfven surface (ie. conical), 
are comparable to stationary solutions obtained with full 
2D, time-dependent nume rical simulations of "cold" jets 
(Krasnopolsky et al. 1999 ). 



Ustyugova et al. (1999), using another MHD code and 
starting with a split-monopole poloidal magnetic field, ob- 
tained a non-conical Alfven surface. However, their local 
jet parameters remain comparable to those obtained in 
this study. Indeed, they have larger mass loads and smaller 
magnetic lever arms, but with a prescribed temperature 
satisfying C+ /nKoTo = {T+ /Tof/'^e > 0.3. Because the 
authors assumed a flat disc (no flaring at all, hence e <C 1), 
this boundary condition requires a coronal temperature 
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larger than the temperature at the disc midplane. As a 
consequence, a significant entropy generation must take 
place in the underlying accretion disc. 



Ouyed & Pudritz (1997) also obtained stationary jets 
from Keplerian discs with full 2D, time-dependent MHD 
simulations. They assumed power-law magnetic fields for 
their boundary conditions at the jet basis and "cold" adi- 
abatic magnetic surfaces. Surprisingly, they obtained jets 
with very large mass loads {k > 0.6) and magnetic lever 
arms very small (A < 3.5). This large mass flux is com- 
patible with a large toroidal magnetic field at the jet basis 
(~ 1.5 times the vertical field). However, the disc vertical 
equilibrium cannot handle such a strong pinching in a non- 
dissipative disc (/ ^ 1, see CF). Besides, we showed in 
Sect. 3.3 that thin discs with T+ < To require A > 4, 
a value smaller than those obtained in these simulations. 
One possible way to reconcile our analytical result with 
their numerical simulation is to realize that the jet basis 
used in all MHD simulations of that kind is not the disc 
surface but some surface above it. In particular, it could 
be located above a region where a coronal heating (allow- 
ing a large mass load k) was followed by a very strong 
cooling (neghgible enthalpy at the "jet basis"). One must 
then be very cautious when identifying boundary condi- 
tions as being the disc surface. Indeed, the real issue of 
jet formation is the mass load provided by the underlying 
accretion disc. 

The effects of entropy generation, theoratically shown 
in this paper, were numerically confirmed by the self- 
similar solutions. For example, the surplus of thermal en- 
ergy, obtained with a significant parameter /, allowed jets 
that would be labelled as slow magnetic rotators {uja < !)■ 
We did not intent to explore any parameter space because 
of the extreme sensibility of the mass load k to the entropy 
source Q (both vertical profile and amplitude, measured 
by /). For example, a very weak entropy generation trig- 
gered inside the disc could give birth to very tenuous, fast 
jets while a huge entropy generation just above the disc 
surface could make dense, slow jets. 

Without exploring the parameter space, we neverthe- 
less obtained a much wider range of MAES parameters 
than under an adiabatic assumption (CF), namely mass 
loads 10~^ < K < 1 (ie. 0.001 < ^ < 0.5), magnetic lever 
arms 1.9 < A < 440 and magnetic rotator parameters 
0.7 < UJA < 2.2. Such huge ejection efficiencies could have 
a strong impact on disc spectral energy distributions. In- 
deed, if the disc is optically thick, the effective tempera- 
ture writes (Ferreira & Pelletier 1995) 



Te// oc r"3 + 4 . (47) 

Thus MAES with large ejection efficiencies would display 
a flat distribution (possibly up to T^fj oc r~^/^). However, 
only a precise treatment of the plasma energetics (includ- 
ing all heating, cooling and transport processes) provides 
the exact value of the disc ejection efficiency ^. This can 
only be done for a given astrophysical context, since all 



Fig. A.l. Cross section of the volume V occupied by the 
disc and the corresponding bounding surfaces (lateral 
and Se, upper S+ and lower S~). The upper and lower 
surfaces correspond to the location of a vanishing radial 
velocity (typically at z+ ~ h{r)). 

these energetic processes (and possibly MHD turbulence 
ctm, Vm and x™) strongly depend on the plasma prop- 
erties. Therefore, any quantitative modelling of a MAES 
must include a complete treatment of the energy equation. 
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Appendix A: Global energy conservation 

The local energy conservation equation can be written as 



V • (/5Wp(y + <i>G + H) + Smhd - u-T) 

pTUp ■ VS - Teff ■ 



(A.l) 



where Smhd is the MHD Poynting flux, T the turbu- 
lent stress tensor, H the enthalpy, S the specific entropy 
and Teff the effective Ohmic and viscous heating (see 
Sect. 3.1). The second law of thermodynamics (p3). 



pTUp ■ VS = Teff + {Tturb 
— V • Srad 



A 



turh 



) + Fe 



(A.2) 



is then incorporated into Eq. (AJ). To get the global en- 
ergy conservation, we integrate this equat ion on the whole 
volum e V occupied by the disc (Ferreira 1994, Henriksen 
1996 ). We thus define S+ and as the disc surfaces 
and SjfSe) the lateral surfaces at r = Vi (r = Vf,, see 
Fig. ( [A.l| )). The disc surface is here precisely defined at 
the locus where the radial velocity vanishes, namely where 
Wr(r, z+ = x^h{r)) = 0. 



After integration, we get 



Pace + Pext + Pt 



turh 



2R 



jet 



(A.3) 



where the accretion power Pace is the power liberated by 
the accretion flow (the difference between what comes 
in at Te and goes out at r^), P^xt = ^yTextdV is the 



power provided by an external source, Pt 



turh 



turh ' 



l^turb)dV is a source of turbulent energy, Pjet is the power 
carried by the outflowing matter and Prad = J^± Srad ■ dS 
is the luminosity emitted at the disc surfaces. 
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The available turbulent power Pturb comes from the 
small scale turbulent motions, but is carried along by the 
laminar accretion flow. Therefore, there is also a flux of 
specific turbulent energy Sturb crossing the inner and outer 
sides of the disc, which implies 



Pturh — 



/ pUpSturb ■ dS pUpSturb ' dS (A. 4) 



We redefine the accretion power Pace by including Pturb 
inside. The total accretion power is then 



ace 



J [pUp£ + SmHD + Srad ~ U ■ J] ■ dS 

/ [pUpE + Smhd+ Srad-u-J] - dS (A.5) 



where the expression of the MHD Poynting flux in the 
resistive regime is given in Ferreira & Pelletier ( |l993aD 
and 



E ^ — +<^>G + H + Eturb 



(A.6) 



is the advected specific energy of the flow (laminar and 
turbulent). The jet power is given by 

Pjet = / [Smhd + pUpE] ■ dS = pUpE{a) ■ dS 

(A.7) 

where E{a) is the Bernoulli invariant ( p7| ) since, at the 
disc surface, the flow is in ideal MHD regime. Finally, we 
obtain 



Pace ^" Pext ^Pjet ^" ^Prad 



(A. 



which is the global energy conservation equation of a 
MAES. 



Appendix B: Self-similar energy conservation 

In the self-similar framework used here, the global en- 
ergy conservation can be expressed in terms of the fiducial 
power 



^ ^ fB.f) 



Plib = Viib 7^ where rjub (B.f) 

2r, f - ^ 

and local disc parameters. The accretion power is a sum 
of various terms, namely Pace = PJ^cT + K'ie + Pa^"° + 
Plcc + -Pace + Pturb- The mcchanical contribution to Pace 
is 



pmec 
riib 



(B.2) 



5 = flo/Qii (f2^ is the Keplerian rotation rate) and rrig — 
Upo/VLxho being defined in CF. The advection of enthalpy 



by the accreting material decreases the liberated energy 
by a factor 



P, 



th 



^2(1-0- 



7 



(B.3) 



Pub ' "7-1" 

while the MHD Poynting flux gives a positive contribution 



pMHD 1 
acc „ , 



e A 



Pi 



lib 



2 1 + A 



B 



qBo 



tan d+S^e 



(B.4) 



where 9^ is the angle between the magnetic field lines 
and the vertical axis at the disc surface {x — x+). The 
presence of a turbulent viscosity produces a radial angular 
momentum flux through the disc, hence a flux of energy. 
This contribution to Pace is 



pvts 

acc 



Pub 

1 -I- A 2rjiib re 



[ iu-T)-dS (B.5) 



where the second term is evaluated at the inner edge of the 
disc. There, it is usually assumed that the viscous torque 
is equal to zero, T(ri) = 0, leaving only a loss of energy 
at the disc outer edge. 

Finally, the total power carried by the jets is 



2R 



jet 



A 



lib 



1 + A 



qBo 



(I -£x+ tan6'+)<5^ 



27 Tt^^. 



1 To 



(B.6) 



(1+ £2^+2)1/2 y^^^^ 

where u is the total velocity of matter. 
Appendix C: Disc thickness 

In this appendix, we focus on the calculation of the disc 
aspect ratio imposed by the energy transport. To illustrate 
the full range of possible astrophysical contexts, we first 
consider systems where the disc plasma is "cold", dense 
and slow (young stellar objects) then systems where the 
disc plasma is "hot" , tenuous and fast (around compact 
objects). 

It is useful to identify which heating process domi- 
nates. In the case of MAES, we are able to estimate the 
ratio of viscous heating to ohmic heating at the disc mid- 
plane, namely 



rjm J 

and 



47'„,M(1 + A)2e2 



3/i 



(C.I) 



(C.2) 



z=0 



For thin [e <C 1) or even slim disc (e < 0.3), we have 
Vm ~ 1, ~ 1, Xm < 1 and A < . The viscous heating 
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is always dominant for solutions with comparable torques 
(A~l). 

We have the same approach than in the standard 



model (Shakura & Sunyaev 1973) since the viscous heat- 
ing is the dominant heating source (we neglect ohmic heat- 
ing) and since the disc is assumed to be optically thick. 
Although the energy transport processes are unknown, we 
can safely assume that the fraction (1 — /) of heating is 
fully transported to the disc surface where it is radiated 
at the photosphere, 



r d7 



dz . 



(C.3) 



The effective temperature T^ff is related to the central 
temperature via the optical depth of the disc [T^fj ~ 

It). Because we have assumed that the total pressure 
gradient equals the plasma pressure gradient, there is a 
direct link between the central temperature To and the 
"photospheric" height of the disc. 

We thus validate the assumption of thin discs in the 
inner regions and checked the assumption of an optically 
thick disc. 

CI. Young Stellar Objects 

In this context, the opacity of the disc can be fitted by 
a Rosseland grain opacity such as nji = O.IT^^^ cm^.g~^ 
(see Bell & Lin|l994|). We obtain 



^phot 



^phot 



0.07(1 a"'/' 



Ma 



2/9 



lO-^Mg/yr / 
The resulting optical depth is 

pKRdz ~ KRPohphot 

0.925 / Ma 



Mq 



-1/3 

(C.4) 



«io/9(l_/)i/9(i + A) \m-^MQ/yr^ 



\1AU 



(C.5) 



This disc is indeed optically thick if the parameter 
and/or the torque ratio A are small {a^ ~ 10~^ in YSO). 

C.2. Compact objects (ACN, microquasars) 

In the inner region of such systems, the Thomson opacity 
dominates. We thus obtain 

£phot — /-i^Q I, J- J ) 



1/10 



M„. 



1/5 



MqIvt \Mq 



-3/10 



r 
3r„ 



1/20 



(C.6) 



where Vg — IGM^Jc? is the Schwarzchild radius. If we 
take M* = lO^M© and Ma = XM^jyv as fiducial values 
for AGN, we get e^hot ^ 3 10-3(1 _ /)i/io„-i/io^ 
In the case of microquasars, we assume Af, = lOAf© 
and Ma = lO-'^Mo/yr, so we get e^hot ^ 10^^(1 - 

The corresponding optical depth is 



2.4 10** 



a^/'^(l-/)i/5(i + A) \3rg 



-3/5 



Ma 



Me/yr 



M, 

Mr. 



-2/5 



3/5 



(C.7) 



As long as the photospheric disc aspect ratio remains 
smaller than unity, these discs remain optically thick. Any- 
way we must be aware that the photospheric height given 
in this appendix is a lower limit since we have neglected 
othe r he ating sour ces. Indeed we can see, by looking at 
Eq.(^T]) and (jcj), that MAES dominated by the mag- 
netic torque (A > e^^) can provide Ohmic heating larger 
(but of the same order) than viscous heating. We thus 
choose a disc aspect ratio larger than those calculated 
here, namely e = 0.1 (a value at the intersection of thin 
and slim discs). 
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